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ABSTRACT 


The  relationship  between  the  number  of  solutions  to  the 
complementarity  problem,  the  right-hand  constant  vector  q 
and  the  matrix  M  is  explored.  The  main  results  proved  in 
this  work  are  summarized  below. 

The  number  of  solutions  to  the  complementarity  problem 
is  finite  for  all  q  c  Rn  if  and  only  if  all  the  principal 
subdeterminants  of  M  are  nonzero.  The  necessary  and 
sufficient  condition  for  this  solution  to  be  unique  for  each 
q  e  Rn  is  that  all  principal  subdeterminants  of  M  are 
strictly  positive.  When  M  >  0  ,  there  is  at  least  one 
complementary  feasible  solution  for  each  q  e  Rn  if  and 
only  if  all  the  diagonal  elements  of  M  are  strictly 
positive;  and,  in  this  case,  the  number  of  these  solutions 
is  an  odd  number  whenever  q  is  nondegenerate.  > .  If  all 
principal  subdeterminants  of  M  are  nonzero,  then  the 
number  of  complementary  feasible  solutions  has  the  same 
parity  (odd  or  even)  for  all  q  c  Rn  which  are  nondegenerate. 
Also,  if  the  number  of  complementary  feasible  solutions  is  a 
constant  eor  each  q  e  Rn  ,  then  that  constant  is  equal  to 
one  and  M  is  a  P-matrix. 

Most  of  the  proofs  are  based  on  mathematical  induction. 
Counterexamples  are  given  to  show  that  the  theorems  fail  if 
any  of  the  hypotheses  are  not  satisfied. 
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ON  THE  NUMBER  OF  SOLUTIONS  TO  THE  COMPLEMENTARY 
QUADRATIC  PROGRAMMING  PROBLEM 

by 

Katta  G.  Murty 

1.  INTRODUCTION: 

1.1  The  complementary  quadratic  programming  problem  is  that  of  finding 
column  vectors  w  «  (w^)  e  Rn  and  z  *  (z^)  e  Rn  satisfying 

w  -  Mz  +  q 

T  (1) 
w  >  0  ,  z  >  0  ,  w  z  ■  0 

where  M  »  (m^)  is  s  given  n  x  n  square  matrix  and  q  ■  (q^)  is  a 

T  n 

given  n  x  1  column  vector  and  w  denotes  the  transpose  of  w  .  R  Is 
the  n-dimensional  real  Euclidean  space. 

1.2  Because  w  ,  z  are  nonnegative,  the  constraint 

T  n 

w  z  ■  y  w.z.  ■  0  =>  w. z.  ■  0  for  each  i  »  1,  ...»  n  . 

!■!  i  1  11 

Thus  if  one  of  the  variables  in  the  pair  w^  ,  z^  Is  positive,  the  other 

T 

should  be  zero.  Hence  the  constraint  w  z  ■  0  will  be  referred  to  as  the 
complementarity  condition  and  the  problem  Is  sometimes  known  as  the 
complementarity  problem  of  order  n  . 

1.3  Consider  the  quadratic  programming  problem 

T 

Minimize  w  z 
Subject  to  w  -  Mz  ■  q 

w  >  0  ,  z  >  0  . 

If  (1)  has  any  solution  (w;z)  then  that  solution  also  solves  the 
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quadratic  programming  problem.  Conversely  if  the  minimum  value  for  the 
objective  function  in  the  quadratic  programming  problem  is  zero,  then  any 
optimal  solution  to  it  also  solves  (1). 

Thus  solving  (1)  is  equivalent  to  finding  out  whether  the  minimum 
objective  value  in  the  above  quadratic  program  is  zero  or  strictly  positive. 
Hence  the  problem  (1)  is  known  as  the  complementary  quadratic  programing 
problem. 

1.4  Cottle  and  Dantzig  [1]  and  Lemke  [7],  [8]  have  shown  that  all  the 
problems  in  linear  programming,  convex  quadratic  programming  and  also  the 
problem  of  finding  a  Nash  equilibrium  point  of  a  bimatrix  game,  can  be  posed 
in  the  form  of  (1).  For  other  applications  of  (1)  see  Scarf  [13].  Lemke 
and  Howson  [7],  [9]  have  developed  a  simple  algorithm  for  solving  (1) 
which  is  based  on  pivot  steps. 

Lemke  [7],  Cottle  and  Dantzig  [1]  have  shown  that  (1)  has  a  solution  if 
all  the  principal  determinants  of  M  are  positive  or  if  M  is  a  nonnagative 
matrix  with  positive  elements  in  the  principal  diagonal.  Lemke  [7]  has  also 
given  sufficient  conditions  on  M  and  q  under  which  the  number  of  solutions 
to  (1)  is  finite. 

In  this  paper  our  main  Interest  is  to  examine  the  relationship  of  the 
number  of  solutions  to  (1)  to  the  properties  of  the  given  matrices  M  and 
q  .  The  motivation  for  this  problem  was  provided  by  Gale  when  he  asked  me 
to  try  and  prove  or  construct  counterexamples  to  the  following  conjectures: 

(a)  M  is  a  P-matrix  if  and  only  if  the  complementarity  problem 
has  a  unique  solution  for  each  q  e'Rn  . 

(b)  M  >  0  ,  is  a  Q-matrix  if  and  only  if  m^  >  0  for  all 


i  -  1 


»  •  •  •  9 


n 


(c)  If  M  is  a  Q-matrix,  the  complementarity  problem  has  an  odd 


number  of  solutions  whenever  q  is  nondegenerate  with  respect 
to  M  . 

The  result  of  the  investigation  is  the  present  work. 
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2.  NOTATION  AND  PRELIMINARIES: 

T 

2.1  If  A  is  any  matrix,  A  denotes  Its  transpose.  A.  denotes  the 
i-th  row  vector  of  A  and  A  .  denotes  the  j-th  column  vector  of  A  .  I 

•j 

denotes  the  unit  matrix. 

2.2  A  square  matrix  M  la  called  a  P  -matrix  If  all  Its  principal  sub- 

determinants  are  strictly  positive.  The  square  matrix  M  Is  called 

nondegenerate  If  every  matrix  A  obtained  by  taking  A  ^  to  be  either 

M  .  or  I  .  for  each  J  -  1,  ...,  n  Is  nonsingular.  An  equivalent 
•J  ‘j 

definition  is  that  M  Is  a  nondegenerate  matrix  If  and  only  if  all  its 
principal  subdeterminanta  ere  nonzero.  M  is  called  a  Q-matrix  if  problem  (1) 
has  a  solution  for  all  q  c  Rn  . 

2.3  Lot  A  be  any  finite  set  of  column  vectors  In  Rn  .  The  convex 
cone  generated  by  the  column  vectors  In  A  Is  denoted  by  Pos  {A}  .  Thus 
x  c  Pos  (A)  if  and  only  If  x  can  be  expressed  as  a  nonnegative  linear 
combination  of  the  column  vectors  in  A  . 

2n 

2.4  Suppose  L(q)C  R  is  the  linear  manifold  determined  by  the 
linear  equality  constraints 


w  -  Mz  -  q 

_n  _n 

v  e  R  ,  z  e  R 


(2) 


without  any  nonnegativity  constraints.  The  vector  (~\~)  c  I»(q)  If  and 
only  if  it  satisfies  (2).  Any  (~\~)  e  L(q)  is  called  a  solution.  For 
convenience  vs  will  write  down  the  vector  (~\~)  e  R2n  **  (wjz)  . 

2.5  The  convex  polyhedron  K(q)  C  L(q)  is  the  set  of  all  feasible 


solutions  (wjz)  which  satisfy 
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w  -  Mz  =  q 
v  >  0  ,  z  >  0 


2.6  A  basic  feasible  solution  is  a  feasible  solution  (w;z)  e  K(q) 
such  that  the  column  vectors  in  (3)  of  the  variables  w,  and  z  which 

J  J 

are  strictly  positive,  are  linearly  independent.  Every  basic  feasible 
solution  is  an  extreme  point  of  the  convex  polyhedron  K(q)  and  vice  versa. 

2.7  A  complementary  feasible  solution  is  a  feasible  solution 
(w{z)  e  K(q)  which  satisfies  the  complementarity  condition  w  z  -  0  . 

A  complementary  feasible  solution  is  a  solution  to  (1)  and  vice  versa. 

2.8  For  each  i  ■  1,  ...,  n  the  variables  w^  ,  z^  constitute  a 

complementary  pair  and  each  of  the  variables  in  the  pair  is  the  complement 

of  the  other.  In  the  system  (1)  the  column  vector  I  is  associated  with 

•J 

the  variable,  w  and  -M  .  is  associated  with  z.  .  Thus  the  pair 
JO  J 

(I  j  ,  -H  j)  are  the  j-th  complementary  pair  of  column  vectors  in  (1). 

2.9  A  complementary  set  of  column  vectors  is  a  set  of  column  vectors 

{A  .  ,  J  ■  1 . n}  such  that  A  .  is  either  I  .  or  -M  .  for  each 

O  O  O  O 

J  ■  1,  . ..,  n  .  Thus  any  set  of  column  vectors  containing  exactly  one  vector 
from  each  complementary  pair  of  vectors  is  a  complementary  set  of  column  vectors. 
The  corresponding  set  of  variables  is  called  a  complementary  set  of  variables. 
Hence  there  are  2n  complementary  sets  of  column  vectors. 

2.10  Each  solution  to  (1)  represents  q  as  a  nonnegative  linear 
combination  of  some  complementary  set  of  column  vectors. 

Conversely  if  {A  .}  is  a  complementary  set  of  column  vectors  and  if 

.  j 

q  e  Pos  (A  ,  ,  j  -  ls  ....  n) 

•J 

n 

i.e.,  q  ■  J  8.A  .  where  Bj  >  0  for  each  j 
M  j“l  J  *J  j  -  J 
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then  a  solution  to  (1)  Is  obtained  by  setting  the  variables  associated 

with  the  column  A  .  equal  to  for  j  ■  1,  . ..,  n  respectively  and 

•J  J 

all  the  other  variables  in  (w;z)  not  in  this  complementary  set  equal  to  zero. 

The  pos  cone  generated  by  any  complementary  set  of  column  vectors 
la  known  as  a  complementary  cone.  Thus  there  are  2n  complementary  cones 
and  the  union  of  all  these  cones  is  the  set  of  all  q  for  which  (1)  has 
a  aolutlon. 

2.11  Any  set  of  variables  {y^  y^.y^ . yQ}  where  yf 

is  either  wf  or  z^,  for  each  r  ,  is  known  as  a  eubcomplementary  set 
of  variables.  The  column  vectors  associated  with  a  subcomplementary  set 
of  variables  constitute  a  subcomplementary  set  of  column  vectors.  The 
complementary  pair  of  variables  (w^z^  is  the  left  out  complementary 

pair  of  variables  in  the  subcomplement  ary  set  {y^ . yi-l,yi+l’  •••»  yn^  * 

2.12  An  almost  complementary  feasible  solution  Is  a  feasible  solution 
(w;z)  e  K(q)  such  that 


w  z  ■  w  z 
11 


for  some  1 


l.e.,  WjZj  ■  0  for  all  J  j*  1  ,  for  some  1  . 

2.13  The  set  (^(q)  Is  the  almost  complementary  set  defined  by 

T 

C1(q)  -  ( (w;z) :  (w;z)  c  K(q)  ,  w  z  ■  w^ 

1...,  Wjtj  -0  for  3^1). 


where  1  Is  any  Integer  from  1  to  n 


2.14 


2.15 


If  x  e  Rn  ,  x  ^  0  then  the  ray  generated  by  x  is 

Pos  {x}  =  {y:  y  =  Xx  for  some  X  >  0}  . 

If  x^"  ,  x^  e  Rn  ,  x1  ^  0  then  the  set 
2  1 

{y:  y  *  x  +  Xx  for  some  X  >  0} 


2  1 
is  the  half -line  through  x  parallel  to  the  ray  generated  by  x 

2.16  The  column  vector  q  is  said  to  be  nondegenerate  with  respect 
to  M  if  and  only  if  for  all  (w;z)  e  L(q)  ,  at  most  n  of  the  2n 
variables  {w^  ,  z^}  are  zero.  Equivalently  q  is  nondegenerate  with  respect 
to  M  if  it  does  not  lie  in  any  subspace  generated  by  (n-1)  or  less 
column  vectors  of  (I  I  -M)  .  Thus  the  set  of  all  q  which  are  not 
nondegenerate  with  respect  to  M  belong  to  a  finite  number  of  subspaces 


of  R 


n 


11  2  2 

2.17  Two  basic  feasible  solutions  (w  ;z  )  and  (w  ;z  )  are  said 

to  be  adjacent  extreme  points  of  K(q)  if  every  convex  combination  of 
11  2  2 

(w  jz  )  and  (w  ;z  )  has  a  unique  representation  as  a  convex  combination 
of  extreme  points  of  K(q)  .  The  line  segments  joining  any  pair  of  adjacent 
extreme  points  of  K(q)  is  called  an  edge  of  K(q)  . 

2.18  If  K(q)  is  nonempty  and  unbounded,  any  basic  feasible  solution  of 


w  - 

n 

+  l 

i-1 


Mz  - 

Zi  " 


0 

1 


w,z  >  0 

m 


is  known  as  an  extreme  homogeneous  solution  of  (3).  Any  half-line  through 
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a  basic  feasible  solution  in  K(q)  parallel  through  the  ray  generated  by 
an  extreme  homogeneous  solution  of  (3),  lies  in  K(q)  .  Such  a  half-line 
Is  called  an  unbounded  edge  (or  extreme  half-line )  of  K(q)  if  every  point 
on  the  half  line  has  a  unique  representation  as  the  sum  of  a  convex  combin¬ 
ation  of  basic  feasible  solutions  of  K(q)  and  a  nonnegative  linear 
combination  of  extreme  homogeneous  solutions  of  (3) . 

2.19  Consider  the  set  of  equality  constraints  (2)  again 

w  ■  Mz  +  q  .  (2) 

The  i-th  constraint  in  this  system  is 

wi  "  Mi.z  +  qi  *  (2D 

A  principal  pivot  in  the  position  (i,i)  in  (2)  consists  of  the  following 

steps: 

(i)  Solve  equation  (2i)  for  the  variable  zi  in  terms  of 

s 

•••*  •••»  zn  an<*  replace  the  i-th  equation 

in  (2)  by  this  equation  expressing  z^  in  terms  of 

*1»  *i-l»wi,*i+l’  zn  * 

(11)  Substitute  the  expression  obtained  for  z^  in  (i)  in  each 
of  the  other  equations  in  (2). 

Thus  a  principal  pivot  in  position  (i,i)  in  (2)  can  only  be  performed  if 
°11  **  ®  ’  The  result  of  this  principal  pivot  is  to  exchange  the  variables 
(v  ,t  >  and  we  get  a  transformed  system  of  equations  which  has  the  same 
form  as  (2),  but  the  left-hand  set  of  variables  in  it  differ  from  the  left- 
hand  set  in  (2)  in  one  component  (the  1-th).  However,  the  set  of  the 
complimentary  pairs  of  variables  remains  unchanged  as  a  result  of  a  principal 
pivot. 
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2.20  If  a  series  of  principal  pivots  are  performed  on  the  system 
(2),  then  it  will  be  transformed  into  the  system 

u  «=  Mv  +  q  (4) 


where  each  pair  (u^,v^)  *s  a  Permutation  of  the  complementary  pair  of 
variables  (w^,z^)  .  A  complementary  feasible  solution  to  (4)  is  a  solution 
to  the  system 


u  *  Mv  +  q 

T 

u>0,v>0,u  v  =  0  . 

e  '  r  * 


(5) 


2.21  We  notice  that  there  is  a  one  to  one  correspondence  between 
solutions  to  (1)  and  solutions  to  (5).  For  example,  suppose  (5)  is 
obtained  from  (1)  by  making  only  one  principal  pivot  in  which  w^  ,  z^ 
are  exchanged,  say.  Then 

A  A  AAA  A 

w  i  z  solves  (1)  <=>  u  *  (z^w^,  ...»  wn>  , 

AAA  A 

v  -  (w,  ,z_,  z  )  solves  (5)  . 

i  i  n 

In  general  since  u  ,  v  in  (5)  are  such  that  *s  a  permjtation 

of  the  variables  (w^,z^)  ,  we  can  construct  a  solution  (ujv)  to  (5) 
corresponding  to  each  solution  (wjz)  to  (1)  by  taking  the  same  permutation, 
and  vice  versa. 

Thus  the  number  of  solutions  to  (1)  is  invariant  under  principal  pivots. 

2.22  Let  N  be  a  principal  submatrix  of  M  of  order  s  ,  obtained 
by  striking  off  from  M  all  the  rows  excepting  the  i.. ,  ...,  i  -th  rows 

X  s 

and  all  but  the  i, ,  ...,  i  -th  columns.  Let 
1  s 


I  Hmmnm 


myr. 


10 


and  Q.  *  /q.  ,  ....  q  Y*  .  Then 

\  1  M 

u  -  N5  +  Q. 

T 

uj0,5>0,iii{»0 


(6) 


Is  known  as  a  principal  subproblem  of  (1)  in  the  variables  (ujO  . 

A  A 

2.23  Suppose  (w;z)  is  a  complementary  feasible  solution  to  (1) 
such  that 


z.  »  0  for  all  i  I*  i.  or  i_,  ....  or  i  . 
l  v  i  s 

Let  a)  -  wi  jT  and  £  -  ^  ,  . . . ,  j  .  Then  from  the 

A  A 

definition  of  the  principal  subproblem  (6)  we  see  that  (w;£)  solves  (6). 

2.2A  If  r  is  any  integer,  its  parity  is  said  to  be  odd  if  r  is 
an  odd  integer  or  even  if  r  is  an  even  integer.  When  considering  a  set 
of  integers,  it  Is  said  to  be  of  constant  parity  if  all  the  numbers  in  the 
set  have  the  same  parity. 

2.25  A  set  of  cones  in  Rn  whose  union  is  R°  is  said  to  form  a 
partition  of  Rn  if  each  cone  in  the  set  has  a  nonempty  interior  and  the 
intersection  of  the  interiors  of  any  two  cones  in  the  set  is  empty. 
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3.  FINITENESS  OF  THE  NUMBER  OF  COMPLEMENTARY  FEASIBLE  SOLUTIONS: 

3.1  Lemke  [7]  has  shown  that  the  number  of  complementary  feasible 
solutions  is  finite  whenever  q  is  nondegenerate  with  respect  to  M  . 
Here  we  determine  the  necessary  and  sufficient  conditions  under  which  the 
number  of  solutions  to  (1)  is  finite  for  each  q  e  Rn  . 

3.2  Theorem: 

The  number  of  complementary  feasible  solutions  is  finite  for  all 
q  e  R°  if  and  only  if  M  is  nondegenerate. 


Proof : 

Suppose  there  exists  a  q  e  R°  such  that  (1)  has  an  infinite 
number  of  distinct  solutions.  Each  solution  to  (1)  represents  q  as  a 
nonnegative  linear  combination  of  some  compl.mentary  set  of  column  vectors. 
There  are  only  2°  distinct  complementary  sets  of  column  vectors.  Thus  if  (1) 
has  an  infinite  number  of  distinct  solutions,  there  must  exist  a  complementary 
set  of  column  vectors  {A  ^  ,  j  ■  1,  . . . ,  n}  such  that 


n 


l 

j-1 


A-JyJ 


■  q 

>  0 


for  each  j  ■  1, 


n 


(7) 


has  an  infinite  number  of  distinct  solutions.  (7)  is  a  square  system  of  n 

equations  in  n  nonnegative  variables.  If  (7)  has  an  infinite  number  of 

solutions,  then  the  set  of  column  vectors  {A  ,  j  ■  1,  ...,  n}  must  be 

•  j 

linearly  dependent.  Since  {A  .}  is  a  complementary  set  of  column  vectors, 

•  J 

this  implies  by  the  definition  in  2.2  that  M  is  not  nondegenerate. 


To  prove  the  converse,  suppose  M  is  not  nondegenerate. 

Case  1:  Suppose  one  of  the  column  vectors  of  M  ,  say  M  ^  ,  is  zero. 

T 

Then  let  q  -  <0,1,1,  ...,  1)  .  Then  (w;z)  *  (0,1,1,  ...,  l;a,0,0 . 0) 

is  a  complementary  feasible  solution  for  any  a  >  0  .  Thus  there  are  an 

infinite  number  of  distinct  complementary  feasible  solutions  when 

T 

q  ■  (0,1,1,  ....  1)  in  this  case. 

Case  2:  Suppose  ft  0  .  Since  M  is  not  nondegenerate,  there 

exists  a  complementary  set  of  columns,  say  (A  ^  ,  j  ■  1,  ....  n}  which 

T 

is  linearly  dependent.  So  there  exists  a  *  (a^ . a^)  i  0  such  that 


l  A.j°j  “  0  * 

j-1  J  J 


Also  A  ,  is  either  I  ,  or  -M  ,  and  hence  in  this  case  A  .  ft  0  . 

• 1  *1  • 1  *1 

n 

If  £  A  .  -  0  ,  let  q  -  A  .  ft  0  .  Then  every  (w;z)  obtained  by 

J-1  '3 

setting  the  variable  associated  with  A  ^  equal  to  1  +  a  ,  the  variable 

associated  with  A.^  equal  to  a  for  J  ft  1  ,  and  all  other  variables  in 

(w;^  equal  to  zero  is  a  complementary  feasible  solution  for  any  a  >  0  . 

Hence  there  are  an  Infinite  number  of  dlstict  complementary  feasible  solutions 

when  q  -  A  ,  ft  0  in  this  case. 

•1 


If  l  A  .  i  0  ,  let  q  -  l  A  .  Let 
J-1  3  J-1  3 

0  -  Min 

J  such  that  \  °J  / 


0  -  Min  [7s-) 

J  such  that  \  °j  / 

Oj  <  0 

-  4«  if  there  does  not  exist  any  <  0  . 

So  0  >  0  .  Every  (w;z)  obtained  by  setting  the  variable  associated  with 


A  j  equal  to  1  +  Xa^  for  j  -  1,  ..., 


n  and  all  the  other  variables  in 
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(w;z)  equal  to  zero,  is  a  complementary  feasible  solution  for  any 

that  0  <  A  <  0  .  Hence  there  are  an  infinite  number  of  distinct 

n 


feasible  solutions  when  q  -  £  A  ^  0  in  this  case 


J-l 


Hence  if  M  is  not  nondegenerate  there  exists  a  q  ^  0  for 
has  an  infinite  number  of  distinct  solutions.  This  completes  the 
Theorem  3.2. 


\  such 
complementary 


which  (1) 
proof  of 


3. 3  Corollary: 

If  M  is  not  nondegenerate  there  exists  a  q  ^  0  for  which  there  are 
infinite  number  of  distinct  complementary  feasible  solutions. 


4.  UNIQUENESS  OF  THE  COMPLEMENTARY  FEASIBLE  SOLUTION; 

4.1  We  will  now  examine  the  question  of  when  (1)  has  a  unique  solution 

for  each  q  c  Rn  while  M  is  fixed.  Lemke  [7],  Cottle  and  Dantzig  [1] 

have  shown  that  (1)  has  a  solution  if  M  is  a  P-matrix.  They  have  also 

shown  that  the  solution  is  unique  if  M  is  positive  definite.  We  will 
extend  these  results. 

4.2  Theorem: 

If  M  is  a  P-oatrix  then  (1)  has  a  unique  solution  for  each 
_n 

q  e  R  . 

Proof : 

Proof  is  by  induction  on  n  .  If  n  ■  1  ,  then  M  ■  (m^)  and  M 
is  a  P-matrix  =i>  >  0  .  In  this  case  q  -  (q^)  .  If  >  0  then 

w  »  (w^  »  (q^  and  z  ■  (z^  ■  (0)  is  the  only  solution  to  (1) .  If 

q^  <  0  ,  then  w  -  (w^)  ■  0  and  z  ■  (z^)  *  (-q^  is  the  only  solution 

to  (1).  So  the  theorem  is  verified  for  n  *  1. 

Suppose  the  theorem  is  true  for  all  complementarity  problems  of  order 
1,2,  ...,  n-1  .  We  will  now  prove  that  it  also  holds  for  problems  of 
order  n  . 

In  Theorem  6  of  [1]  Cottle  and  Dantzig  have  proved  that  (1)  has  at 
least  one  solution  when  M  is  a  P-matrix.  Parsons  has  proved  the  same 
by  using  induction  on  n  (Theorem  6.1  of  [11]).  Thus  we  need  only  prove 

A  A 

the  uniqueness  of  the  solution.  Suppose  a  solution  to  (1)  is  (w;z) 


where 
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z  >  0,  . . . ,  z  >0 
J1  Jr 

and  (8) 

* 

Zj  ■  0  for  all  j  t  or,  ...»  or  Jr  • 


Make  a  principal  pivot  exchanging  z  and  w.  which  Is  possible  because 

h  J1 

m.  .  >  0  by  the  P-  property  of  M  .  The  result  of  this  principal  pivot 

J1J1 

is  another  system  of  the  same  form  as  (1),  in  which  the  matrix  M*  is 

again  a  P- matrix  by  Tucker's  theorem  [15]. 

In  the  new  system  make  a  principal  pivot  exchanging  z  and  w 

j2  J2 

which  is  possible  because  m'  .  >  0  by  the  P  -  property  of  M'  . 

J2J2 

Continue  making  principal  pivots  until  all  the  original  z  ,  z 

J1  Jr 

are  exchanged  with  w  ,  ....  w  respectively.  This  is  possible  because 

J1  Jr 

after  each  principal  pivot  we  get  a  new  system  in  which  the  matrix  is  again 
a  P-  matrix  by  Tucker's  thereom  [19  and  hence  has  strictly  positive 
elements  along  the  principal  diagonal. 

Suppose  the  system  obtained  at  the  end  is 


u  >  0 


u  ■  Mv  +  q 

T 

,  v  >  0  ,  u  v 


0 


(9) 


By  the  manner  in  which  (9)  was  obtained  from  (1)  and  since  (w;z)  is  a 
solution  to  (1)  we  note  using  (8)  and  2.21  that  (9)  has  a  solution  in  which 
v  •  0  .  So  in  that  solution  u  ■  q  and  since  u  must  be  nonnegative, 
q  >  0  .  Thus  by  a  series  of  principal  pivots  we  have  transformed  (1)  into 
(9)  in  which  M  is  a  P  -matrix  (by  Tucker's  theorem  [L5])  and  q  >  0  . 
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Since  M  is  a  P  -  matrix  all  its  principal  submatrices  are  P~ 
matrices.  Hence  by  the  induction  hypothesis  all  principal  subproblems 
(defined  in  2.22)  of  (  9  )  have  unique  solutions.  (9)  has  one  solution 
u  -  q  ,  v  -  0  .  If  (9)  has  another  distinct  solution  (u*;v*)  in  which 
for  some  J  ,  v*  ■  0  ,  then  the  principal  subproblera  of  (9)  in  the 

variables  (u^,  ...,  •••»  un  »  vi>  •••»  vj-i,vj+i’  •••»  vn^ 

has  two  distinct  solutions ,  namely 


*  q^  for  e»ch  i  i  j  and  ■  0  for  each  i  +  J 

and 

u.  -  u*  for  each  i  4  j  and  v1  «  v*  for  each  i  +  j 

and  these  two  are  distinct  because  (u*;v*>  is  distinct  from  (q;0)  . 
This  contradicts  the  induction  hypothesis. 

Hence  if  (  9)  has  another  solution  (u*;v*)  distinct  from  (v;0) 
then  we  must  have  v*  >  0  .  By  the  complementarity  constraint 
v*  >  0  ■>  u*  ■  0  .  Therefore 


Mv*  •  -q  <  0 

^  at 

v*  >  0  . 


(10) 


This  is  a  contradiction  because  by  the  theorem  of  Gale  and  Nikaido 
(Theorem  1  of  (6l)  since  M  is  a  P- matrix. 


Mv  <  0 
v  >  0 


*>  v  ■  0  . 


Hence  (9)  has  a  unique  solution.  By  2.21  we  conclude  that  the  solution  to 
(1)  is  also  unique. 

Hence  by  induction,  Theorem  4.2  is  true  for  all  n  . 
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It  is  interesting  to  note  that  the  converse  of  Theorem  4.2  is  also 
true.  This  is  proved  next. 

4.3  Theorem: 

If  (1)  has  at  most  one  solution  for  each  q  e  Rn  ,  then  M  is  a 
P  -matrix. 


Proof : 

Proof  is  by  induction  on  n  .  Suppose  n  ■  1  .  Then  M  -  (m^) 
and  q  -  (q^  . 

If  ■  0  ,  then  if  q  -  (q^)  ■  0  ,  (1)  has  several  solutions, 
namely 

w  -  (w^  -  0  ,  z  ■  (z^  for  any  z^  >  0  . 

If  Ojj  <  0  ,  then  if  q  »  (q^)  >  0  ,  (1)  has  two  distinct  solutions, 
namely, 


w  -  (qx)  ,  z  -  0 


Thus  if  n  •  1  ,  by  the  hypothesis  of  Theorem  2  m^  £  0  .  So  we  must 
have  m^  >  0  and  hence  M  ■  (m^)  must  be  a  P-  matrix,  which  proves  the 
theorem  for  the  case  n  ■  1  .  Suppose  the  theorem  holds  for  all  problems  of 
order  (n  -  1)  or  less.  Ve  will  now  show  that  it  also  holds  for  problems 


of  order  n  . 
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If  (1)  has  at  most  one  solution  for  each  q  e  Rn  ,  then  every  principal 

subproblem  of  (1),  has  at  most  one  solution  for  any  of  its  right-hand  side 

constant  vector  Q,  .  To  prove  this  consider  the  subprobleir  in  the  variables 
T  T 

W  -  (w_ ,  ....  w  )  ,  5  -  (Z,.  ...»*)  •  If  there  exists  a 

z  it  z  n 

Q.m  Cq2»  •••»  qR)T  for  which  this  principal  subproblem  has  two  distinct 
solutions,  namely,  (ui;0  and  (ui;0  ,  choose  to  satisfy 

ql  >  MaX  Smljl  ’  ^l2  Cjmlj'}  * 


and  let 


n  _ 


'l  '  ql  +  3l2  Ej"lj  1  *1 


-  0 


n  a 


0  . 


"1  '  "l  +  }l2  EJBU  1  'l 

If  q  -  ,  then  (1)  has  two  distinct  solutions,  namely 

5  ‘  (it)  ’ 5 '  ('!”) 


and 


•  (•?•)  •  ■  ■  (f) 


contradicting  the  hypothesis. 

A  similar  proof  holds  for  all  principal  subproblems  of  (1)  of  order 
n  -  3,  ,  and  by  repeating  the  argument  we  see  that  all  principal  subproblems 
of  (1)  have  at  most  one  solution  for  all  right-hand  side  constant  vectors. 
Hence  by  the  induction  hypothesis  all  principal  submatrices  of  M  of 
order  (n  -  1)  or  less  are  P-  matrices.  In  particular  m^  >  0  for  all 


The  matrix  M  itself  should  be  nonsingular,  as  otherwise  it  will  not 
be  nondegenerate,  in  which  case  there  exists  a  q  for  which  (1)  has  an 
infinite  number  of  solutions  by  Theorem  3.2.  So  M  *  exists.  Left 
multiplying  by  M  ^  problem  (1)  can  be  written  as 

-1 

z  -  M  w  »  Q 

T  (11) 

z>0,w>0,zw«0 

where  Q  ■  -M  . 

If  (1)  has  at  most  one  solution  for  each  q  e  Rn  ,  then  (11)  has  at 
most  one  solution  for  each  Q  e  Rn  .  Hence  by  the  arguments  used  previously, 
all  principal  subdeterminants  of  M  ^  of  order  (n  -  1)  or  less  are 
strictly  positive.  Let  a  be  the  value  of  the  principal  subdeterminant 
of  M  *  obtained  by  striking  off  the  first  row  and  column  from  M  *  .  Then 


a 


mll _ 

determinant  of  M 


(12) 


But  a  >  0  ,  >  0  .  So  by  (12)  the  determinant  of  M  is  positive. 

Also  every  principal  submatrix  of  M  is  a  P-  matrix.  Hence  M  Itself 
is  a  P -matrix.  Hence  by  induction  Theorem  4.3  holds  for  all  n  . 

4.4  Corollary; 

If  (1)  has  at  most  one  solution  for  each  q  e  Rn  (M  being  fixed) 
then  it  has  exactly  one  solution  for  any  q  e  Rn  . 

4.5  Corollary: 

(1)  has  a  unique  solution  for  each  q  c  Rn  if  and  only  if  M  is 


a  P  -matrix. 


A. 6  Corollary : 


If  (1)  has  a  unique  solution  for  each  q  e  Rn  ,  then  any  principal 
subproblem  of  (1)  has  a  unique  solution  for  each  of  its  right-hand  side 
constant  vector. 

A. 7  Note;  It  is  not  possible  to  generalize  Corollary  A. 6  by  dropping 
the  condition  of  "uniqueness"  of  the  solution  to  (1).  As  an  example,  let 
M  ■  (  \  .  Then  the  problem  is  to  solve 


.ixy-t) 


w>0,z>0,wz 


(13) 


In  the  diagram  below,  we  indicate  each  complementary  cone  by  drawing 
a  dotted  line  segment  running  across  its  generators. 


We  see  that  the  union  of  the  complementary  cones  is  the  whole  space,  R 

2 

So  by  2.10,  (13)  has  a  solution  for  all  q  e  R  .  So  M  is  a  Q-matrix. 
Now  let  us  examine  the  principal  subproblem  of  (13)  in  ^2  *  z2  ' 


It  is 
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v2  -  (-l)z2  =  q2 
w2  >  0  ,  z2  >  0  ,  «2z2  -  0 


(14) 


(14)  has  no  solution  if  q2  <  0  .  Also  we  verify  that  M  violates  the 
hypothesis  of  Corollary  4.6  because  it  has  three  distinct  Solutions  for  every 
q  in  the  positive  orthant. 

4.8  Thus  the  fact  that  M  is  a  Q-matrix  does  not  imply  that  its 
principal  submatrices  are  Q-matrices. 

4.9  Note:  After  Theorems  4.2,  4.3  were  conjectured  and  proved  as 
described  here,  a  theorem  by  Smelson,  Thrall  and  Wesler  [12]  on  the  partition 
of  Rn  by  convex  cones  which  is  equivalent  to  Corollary  4.5  has  come  to 

our  notice.  Their  proof  is  geometric  and  not  based  on  the  mathematical 
programming  approach. 

The  property  of  "uniqueness"  of  the  solution  to  (1)  also  affects  the 
nature  of  the  solution.  This  is  discussed  below. 


4.10  Theorem: 

Suppose  M  has  the  property  that  (1)  has  a  unique  solution  for  each 
q  e  Rn  .  Keep  q2>  . . . ,  qQ  fixed  but  let  q^  vary.  Let  z^(q^)  be  the 
value  of  z^  in  the  solution  to  (1)  as  a  function  of  q^  .  Then  z^(q^) 
is  monotonic  decreasing  in  q^  and  it  is  strictly  monotonic  decreasing  in 
the  region  in  which  it  is  positive. 

Proof : 

T 

Proof  by  contradiction.  Let  Q.  ■  (q2»  ...,  qfl)  which  is  held  fixed. 

A  A 

Pick  any  value  for  q^  and  let  8  >  0  be  arbitrary.  Let  (w;z)  be  the 

,  (w;z)  be  the  solution  to  (1)  when 


solution  to  (1)  when  q 
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(--4--)  • 


Then  z^q^)  “  z\  and  zl^l  +  8)  “  •  Let 


A  A  A  a  A  A  J 

111  *  (w^  y  •  •  •  |  W  )  I  ^  c  (Za  J  •  •  •  f  z  ) 

Z  n  z  n 

-  -  •£  -  t 

u  ■  (w0,  ...,  w  )  ,  £  =  (z,,  •••»  z  ) 


and 


m#1  -  (®21 ’  •••»  mnl^  * 


If  z^  >  0  we  wish  to  show  that  z\  <  z\  '  Suppose  not,  then 
z^  >  z^  >  0  .  By  complementarity  w^  *  w^  *  0  .  Then,  if 


Ai  + 6  +  vA  _ 

\  «+".i'i  /  ’ 


(1)  has  two  solutions,  namely 


.  -  (-1-)  ,  3  -  M-) 


and 


contradicting  the  uniqueness  of  the  solution  to  (1)  for  each  q  e  Rn  .  So  if 
Zi(qi>  >  0  then  z1(q1  +  B)  <  z^q^  for  any  B  >  0  . 

It  remains  to  be  shown  that  if  z^(q^)  ■  0  then  z^(q^  +  B)  ■  0  for 

A  “* 

all  8  >  0  .  Suppose  not,  then  z^q^)  “  m  0  and  z^4^  +  8)  -  >  0  . 

■  »  (1)  has  two  solutions,  namely 


Then  if  q 


w  ;  z  and 


which  is  again  a  contradiction. 
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A.  11  We  now  show  that  if  M  is  a  Q-matrix  and  (1)  has  a  unique 
solution  when  q  is  any  element  of  the  set 

{I  ,, I  I  ;  -M  . . -M  }  ,  then  M  is  a  P-matrix. 

t i  • t  *n  *1  •  n 

4.12  Theorem: 

Let 


union  of  all  complementary  cones  which  contain 
generator 

union  of  all  complementary  cones  which  contain 
generator. 


as  a 


as  a 


If  l  j  i  [Zj]  and  -M  ^  i  [w^  ]  for  each  j  -  1  to  n  and  M  is  a 
Q-matrix,  then  M  is  a  P-matrix. 

Proof: 

4.13  Let  N  be  the  principal  submatrix  of  M  of  order  (n  -  1)  , 
obtained  by  striking  off  the  first  row  and  column  of  M  .  We  will  now  show 
that  N  is  also  a  Q-matrix. 

T  T 

Suppose  not.  Let  m  ■  (w£,  ...»  wr)  and  £  ■  (*2 . *n)  •  Consider 

the  principal  subproblem  in  (u>,0  ,  which  is  to  solve 

u  -  NC  -  £ 

T  (15) 

w>0,£>0,u)£*0. 

If  N  is  not  a  Q-matrix  there  exists  a  £  e  Rn  *  such  that  when  Q  ■  , 

(I5)  has  no  solution. 

Let  q  -  ^---j  .  If  (1)  has  a  solution  (w,z)  when  q  -  q  ,  then 
I  >  0  in  it,  as  otherwise  (ui,£)  would  be  a  solution  to  (15).  Thus  every 
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point  on  the  line 

f  = 5  •  (-{)  • a  realj 

corresponds  to  only  complementary  feasible  solutions  in  which  Zj  >  0  .  Since 
there  are  only  a  finite  number  of  complementary  cones  and  each  one  is  convex, 
there  must  exist  an  aQ  such  that  the  half-line 


lies  entirely  in  a  complementary  cone.  By  the  above  argument,  in  every 
complementary  feasible  solution  corresponding  to  any  point  on  this  half-line 
we  must  have  zl  >  0  .  This  implies  that  this  half-line  lies  in  a  complement¬ 
ary  cone  for  which  -M^  is  a  generator.  This  implies  that  PosU.^  also 
lies  in  the  same  complementary  cone,  i.e.,  I#1  e  [z^]  >  which  is  a  contradiction 

So  N  must  be  a  Q-oatrix.  By  a  similar  argument  we  conclude  that  all 
principal  submatrices  of  M  of  order  (n  -  1)  must  be  Q-matrices. 

4.14  Let  N  be  the  principal  submatrix  of  M  of  order  (n  -  2) 
obtained  by  striking  off  the  first  two  rows  and  columns  from  M  .  We  will 

now  show  that  N  must  be  a  Q-matrix  also. 

Suppose  not.  Then  there  exists  e  5  c  Rn'2  such  that  the  subproblem  in 
(w  „n  ;  z3,  ...,  zn)  has  no  complementary  feasible  solution  when  its 

right  hand  constant  vector  is  Q  . 


Then  for  any  A,y  real,  (1)  has  no  solutions  in  which  both  z1  and  z2 
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are  equal  to  zero,  when  q  ■  q  by  2.23. 

4.15  Fix  A  =  A^,  y  *  y^  where  A^  >  0  and  y^  >  0  and  consider  the 

line 


Points  on  this  line  have  complementary  feasible  solutions  in  which  both  z^ 
and  z 2  cannot  be  zero  together.  Since  the  number  of  complementary  cones 
is  finite  and  each  is  convex,  there  must  exist  an  a  such  that  the  entire 
half-line 


is  in  a  complementary  cone.  Suppose  this  half-line  is  in  the  complementary 

cone  Pos{I  ,  ,-M  A  A  }  .  Then 

•l  •  j  *n 

A.I  ,  +  y.I  _  e  Pos{I  -,-M  olA  . . A  }  . 

l*x  i  *z  • i  • j  *n 


Suppose 

V.X  +  “Ir.2  ■  V.l  +  +  l3  V-J 

where  a,,ot0,  . i  >  0  . 

12  n  ■ 

If  >  X^  ,  then  if  we  put  (A,y)  ■  (0,y^)  ,  q  of  (16)  will  have  a 
complementary  feasible  solution  in  which  both  z^  -  ■  0  ,  which  is  a 

contradiction  to  4.14. 

If  o,  <  X,  ,  then  (X,  -  a,)I  ,  +  y.I  »  lies  in  t  ■=  intersection  of 
11  1  1  *1  1  ’2 

Pos{-M  0,A  A  }  with  Pos{I  .  ,1  ,}  .  We  note  that  Pos{I  .,1  -} 
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cannot  entirely  lie  in  Posf-M  ^A. -j>  •••»  A  because  then  I  ^  e  [z^] 
and  e  [*2^  contradicting  the  hypothesis.  So  Pos{I  ^,1^  and 

Pos{-M  -,A  ,,...,  A  }  intersect  in  a  half-line  and  when  (A,p)  -»  0  unless 
A  ,p  are  such  that 


A  > 


Xl- 


we  have 


AIt  ^  wl. 2  ^  Pos { I ^ ^ M ^ j »A^ 3 ,  •  ••»  At^}  * 

Similarly  we  see  that  for  (A,p)  >  0  if  AI  ^  +  pl^  is  contained  in 
a  complementary  cone  Pos{-M  ^,1  2*B.3»  •••»  B.n^  »  t*ien  +  yl  2  must 

lie  on  some  half-line  in  Pos{I  ,  ,1  ,}  . 

•  1  •  L 

Hence  when  (A.,y.)  >  0  ,  unless  (A. I  ,  +  p,I  S)  lies  in  the  union  of 

a  finite  number  of  half-lines  in  PosU^.I^}  ,  the  half-line  in  (17)  can 

only  be  contained  in  a  complementary  cone  for  which  both  -M  ^  and  -M  ^ 

are  generators.  This  implies  that  all  the  points  in  Pos {I  ^,1^}  excepting 

those  lying  on  a  finite  number  of  half-lines  belong  to  the  union  of  all 

complementary  cones  containing  both  -M  ^  and  -M  ^  as  generators.  But 

this  union  is  a  closed  cone  and  if  it  contains  all  points  of  Pos{I  ^,1^} 

excepting  those  lying  on  a  finite  number  of  half-lines,  then  it  contains 

all  of  Pos  {I  2^1.2^  •  This  implies  that  I  ^  lies  in  some  complementary 

cone  which  has  both  -M  ,  and  -M  «  as  generators,  which  is  a  contradicti:- 

•  1  ♦  i 

to  the  hypothesis. 

So  N  oust  be  a  Q-matrlx.  By  a  similar  argument  we  can  show  that 
every  principal  submatrix  of  M  is  a  Q-matrix.  Hence  all  the  elements  in  the 
principal  diagonal  of  M  must  be  strictly  positive. 


From  the  hypothesis  of  the  theorem  we  see  that  every  matrix  M  ,  obtained 
from  M  by  performing  a  series  of  principal  pivots  (as  in  2.20)  has  the 
property  that  all  its  diagonal  elements  are  strictly  positive. 

By  Tucker's  theorem  [15]  (see  also  Lemma  6.1  in  [11])  this  implies 
that  M  is  a  P-matrix. 

4.16  Note :  It  may  be  possible  to  use  Theorem  4.12  to  develop  an 
efficient  algorithm  for  testing  whether  a  given  real  square  matrix  M  is 


a  P-matrix  or  not. 
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5.  ON  THE  Q-  NATURE  OF  NONNEGATIVE  MATRICES 

5.1  Suppose  the  square  matrix  M  is  nonnegative,  i.e.,  m^  >  0 
for  each  i  and  j  .  This  case  is  of  particular  interest  because  the 
problem  of  finding  a  Nash  equilibrium  point  of  a  bimatrix  game  can  be 
formulated  as  a  problem  of  the  form  (1)  in  which  M  >  0  .  See  [1].  It 
is  of  interest  to  know  when  such  a  matrix  is  a  Q-  matrix.  The  following 
theorem  discusses  this  question. 

5. 2  Theorem: 

Let  M  >  0  .  M  is  a  Q- matrix  if  and  only  if  >  0  for  each 

i  ■  1,  . . . ,  n  . 

Proof  of  Sufficiency: 

This  has  been  proved  by  Cottle  and  Dantzig  [1]  in  a  corollary  under 
their  Theorem  5. 

Proof  of  Necessity: 

Proof  by  _uduction.  Suppose  n  ■  1  .  If  M  -  (m^)  »  0  then  for 
q  ■  (q^)  <  0  ,  (1)  has  no  solution.  So  the  theorem  is  verified  for  n  ■  1 

Suppose  the  conditions  of  the  theorem  are  necessary  for  all  problems 
of  order  (n  -  1)  or  less.  We  will  show  that  they  are  then  necessary  for 
n  also.  Let  M  >  0  be  a  square  matrix  of  order  n  ,  with  at  least  one 
nii  e<lua*  to  zero*  Without  any  loss  of  generality  we  can  assume  that 


Consider  the  principal  subproblem  of  (1)  in  the  variables 

T  T 

u  ■  (w. i  • • • ,  w  . )  ,  £  -  (z. ,  •••»  ^  i  )  • 

l  n-i  l  n-x 
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The  matrix  corresponding  to  this  subproblem  is  nonnegative  with  a  zero 
element  (m^  *  0)  in  its  diagonal.  So  by  the  induction  hypothesis,  there 
exists  a  Q  *  (q^,  qn_^)  such  that  this  principal  subproblem  of  (1) 

has  no  solution  when  its  right-hand  side  constant  vector  is  Q.  .  Let 


q  >0  and  q  * 


(t> ' 


When  q  ■  q  ,  (1)  has  no  solution  because 


(1)  If  (1)  has  a  solution  at  all,  say  (w  ;  z)  then  z^  +  0  ,  since 

z^  -  0  Implies  that  (u  ;  O  solve  the  principal  subproblem 

contradicting  the  manner  in  which  §.  was  obtained. 

(ii)  So  z  >  0  .  But  by  (1)  w  ■  M  z  +  q  >0  because 
n  n  n»  ’n 


M  >0,z>0,q  > 

n*  -  -  mi 


So  if  (1)  has  any  solution  at  all,  both  w^  ,  zq  must  be  positive  in  it 
which  violates  the  complementarity  constraint. 

Then  (1)  has  no  solution  when  q  is  obtained  as  above.  So  the  theorem 
is  true  for  n  .  By  induction,  it  is  true  for  all  n  . 


5.3  Proof  of  Theorem  5.2  by  Example; 


Suppose  m. .  "0  .  Then  if 


-  (-1,1.1 . 1) 


Then  using  the  fact  that  M  >  0  it  is  easily  verified  that  (1)  has  no 


solution  when  q  ■  q  .  This  example  is  given  by  Gale. 


5. A  Corollary: 


If  M  >  0  and  M  is  a  Q-matrix,  then  all  principal  submatrices  of  M 
are  also  Q-matrlces. 


i 
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6.  ON  THE  CONSTANT  PARITY  OF  THE  NUMBER  OF  COMPLEMENTARY 
FEASIBLE  SOruliON,?: 

6.1  We  will  now  examine  how  the  number  of  solutions  to  (1)  varies 
as  q  varies  over  Rn  while  M  is  fixed. 

6. 2  Theorem: 

If  M  is  nondegenerate,  then  the  number  of  complementary  feasible 
solutions  has  the  same  parity  for  all  q  e  Rn  which  are  nondegenerate 
with  respect  to  M  . 

Proof ; 

In  the  proof  of  this  theorem  we  will  use  some  of  the  results  proved 
by  Lemeke  in  [7], 

6.3  Results  [Lemkel: 

If  q  is  nondegenerate  with  respect  to  M 

(i)  Then  (1)  has  a  finite  number  of  solutions  when  q  «  q  . 

(ii)  For  each  i  -  1,  ....  n  ,  the  almost  complementary  set 

C^(q)  is  either  empty  or  is  the  union  of  some  edges  (bounded 
or  unbounded)  of  K(q)  . 

(iii)  The  number  of  unbounded  edges  in  C^(q)  differs  from  the  number 
of  solutions  to  (1)  by  an  even  number. 

6. A  We  will  now  prove  that  if  M  is  nondegenerate  and  q  e  Rn  then 
is  unbounded  on  every  unbounded  edge  of  C^(q)  .  Suppose  F  is  an 
unbounded  edge  of  K(q)  contained  in  C^(q)  .  Let 

F  -  {(w;z):  (wjz)  -  (w1  +  8w2  ;  z1  +  0z2)  ,  0  >  0} 
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11  2  2 
where  (w  ;z  )  Is  a  basic  feasible  solution  and  (w  ;z  )  is  an  extreme 

homogeneous  solution  of  K(q)  . 

1  2 

Along  this  unbounded  edge  F  ,  +  0w^  ,  and  it  remains  bounded  for 

2  2  '  1 
all  6  >  0  only  if  w^  *  0  .  If  w^  ■  0  ,  then  if  we  put  q  -  q  -  w^I  , 

(1)  has  an  infinite  number  of  solutions,  namely 

w  =  w*  -  Wjl  +  9w2  ;  z  «  z^  +  0z2  for  all  0  >  0 

i  •  i  “ 


which  is  a  contradiction  to  the  hypothesis  that  M  is  nondegenerate  by 
Theorem  3.2. 

Thus,  on  every  unbounded  edge  of  C^(q)  ,  w^  is  unbounded. 

6.5  We  will  now  use  the  result  obtained  in  6.4  to  show  that  if 
q  e  Rn  is  nondegenerate  with  respect  to  M  and  a  is  any  real  number  such 
that  q  *  q  +  al#i  is  also  nondegenerate  with  respect  to  M  ,  then  the 

*  A 

number  of  unbounded  edges  in  C^(q)  and  C^(q)  are  the  same. 

Pick  any  unbounded  edge 


F  ■  {(w;z)  :  (w;z)  «  (w"  +  0w2  ;  z1  +  0z2)  ,  0  >  0}  C  C^q)  . 


By  6.4,  w^  >  0  and  w^  is  unbounded  on  this  edge.  Let 


v  -  Max 


0  , 


1  . 

w.^  +  a 


-w. 


Then  it  is  easily  verified  that 

F*  ■  {(w;z)  :  (w;z)  -  (w*  +  0w2  +  ai#i  ;  z*  +  0z2)  ,  6  >  v} 
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is  an  unbounded  f  .  ...  C^(q)  . 

Thus  we  have  shown  that  there  exists  an  unbounded  edge  in 

C^(q)  corresponding  to  each  unbounded  edge  F  in  C^(q)  .  Conversely 

by  treating  q  E  q  +  (-a)  I  .  we  can  establish  a  correspondence  between 

•  1 

unbounded  edges  in  C^(q)  and  those  of  C^(q)  .  This  establishes  a 
1-1  correspondence  between  the  unbounded  edges  in  C^(q)  and  those  in 
C^(q)  .  Hence  both  C^(q)  and  C^(q)  must  have  the  same  number  of 
unbounded  edges. 

6.6  Now  to  continue  the  proof  of  Theorem  6.2,  let  q  and  q  be  any 
two  column  vectors  in  Rn  both  of  which  are  nondegenerate  with  respect 

to  M  . 

By  6.5  and  (i),  (iii)  of  6.3  we  conclude  that  the  parity  of  the 
number  of  solutions  to  (1)  does  not  change  if  we  alter  the  vector  q  one 
component  at  a  time  so  that  it  remains  nondegenerate  with  respect  to  M 
both  before  and  after  the  alteration. 

It  is  always  possible  to  alter  q  ,  by  one  component  at  a  time, 
retaining  the  property  of  being  nondegenerate  with  respect  to  M  throughout, 
until  it  becomes  equal  to  q  . 

Hence  the  number  of  solutions  to  (1)  has  the  same  parity  whether 
q  ■  q  or  q  .  Hence  the  number  of  solutions  to  (1)  has  the  same  parity 
whenever  q  is  nondegenerate  with  respect  to  M  . 

6.7  Note:  The  assumption  that  M  is  nondegenerate  cannot  be  dropped 
from  the  hypothesis  of  Theorem  6.2,  as  can  be  seen  from  the  example  below. 

Let  M  ■  2)  •  Then  (1)  is  to  solve 


0 


0 


W1  ’  W2  ’  Z1  1  Z2  =  U  ’  W1Z1  +  W2Z2  *  ° 


q  -  leads  to  one  solution  to  (1)  and  q  ■  leads  to  two  solutions 

to  (1)  eventhough  both  these  are  nondegenerate  with  respect  to  M  here. 

Here  M  is  not  nondegenerate  because  the  matrix  is 

singular.  The  argument  in  6.4  fails. 

6.8  Note:  The  assumption  that  q  is  nondegenerate  with  respect  to 
M  cannot  be  dropped  from  the  hypothesis  of  Theorem  6.2  as  can  be  seen  from 
the  example  below.  Let  M  -  Q  • 
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When  q  is  nondegenerate  with  respect  to  M  the  number  of  solutions  to 
(1)  is  an  odd  number,  but  when  q  *  -M  2  *  (1)  has  exactly  two  distinct 
solutions. 

6.9  Corollary: 

If  M  is  nondegenerate  and  not  a  Q-  matrix  then  the  number  of 
solutions  to  (1)  is  an  even  number  for  all  q  which  are  nondegenerate 
with  respect  to  M  . 

6.10  Proof : 

By  2.  10,  the  set  of  all  q  for  which  (1)  has  a  solution  is  the 
union  of  the  2°  complementary  cones.  Each  complementary  cone  is  a  closed 
set  in  Rn  and  hence  their  union  (being  a  union  of  a  finite  number  of 
closed  sets)  is  itself  closed.  The  set  of  all  q  for  which  (1)  has  no 
solution  is  the  complement  of  this  union,  and  hence  is  an  open  set.  Because 
M  is  not  a  Q- matrix,  this  open  set  is  nonempty.  So  the  set  of  all  q 
for  which  (1)  has  no  solution  is  a  nonempty  open  cone. 
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6.11  By  2.16,  the  set  of  all  q  which  are  not  nondegenerate  with 
respect  to  M  is  the  union  of  a  finite  number  of  subspaces  of  Rn  ,  each 
of  which  has  dimension  <  (n-1)  .  Hence  the  set  of  all  q  which  are  not 
nondegenerate  with  respect  to  M  has  no  interior. 

6.12  By  6.10  and  6.11  we  conclude  that  there  must  exist  a  q 
nondegenerate  with  respect  to  M  ,  for  which  (1)  has  no  solution,  i.e., 
zero  solutions.  Now  by  applying  Theorem  6.2  we  conclude  that  the  number 

of  solutions  to  (1)  has  the  same  parity  as  zero,  i.e.,  even  parity,  whenever 
q  is  nondegenerate  with  respect  to  M  . 

6.13  Note:  Corollary  6.9  is  not  necessarily  true  if  M  is  not 
nondegenerate  as  seen  from  Example  6.7. 

6.14  Note:  The  converse  of  Corollary  6.9  is  not  necessarily  true 
unless  M  >  0  .  This  is  discussed  in  8.17. 


:.'-***W  *•  I.*.  MM  M 
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7.  ON  PROBLEMS  WITH  A  CONSTANT  NUMBER  OF  COMPLEMENTARY  FEASIBLE  SOLUTIONS 


7.1  Here  we  show  that  if  the  number  of  complementary  feasible  solutions 
is  a  constant  for  all  nonzero  q  e  Rn  then  that  constant  is  equal  to  one  and 
M  is  a  P-matrix. 

7.2  Theorem: 

If  the  number  of  complementary  feasible  solutions  is  a  constant  for  all 
q  c  Rn  ,  q  |l  0  ,  then  M  is  a  P-matrix  and  that  constant  is  equal  to  one. 

Proof : 

7.3  Whatever  M  may  be,  (1)  always  has  at  least  one  solution  for  every 
q  >  0  (the  solution  is  w  ■  q  ;  z  *  0)  .  If  M  is  not  a  Q-matrix,  there 
exists  a  q  /  0  for  which  (1)  has  no  solution  at  all.  Hence,  if  M  is  not 
a  Q-matrix  then  the  number  of  solutions  to  (1)  cannot  be  a  constant  for  all 

q  )>  0  .  So  under  the  hypothesis  of  Theorem  7.2,  M  must  be  a  Q-matrix. 

7. A  The  number  of  complementary  feasible  solutions  is  finite  whenever 
q  Is  nondegenerate  with  respect  to  M  .  By  Corollary  3.3,  if  M  is  not 
nondegenerate,  there  exists  a  q  f  0  for  which  (1)  has  an  infinite  number  of 
solutions. 

Since  the  number  of  solutions  to  (1)  is  a  constant  for  any  q  +  0  ,  M 
must  therefore  be  a  nondegenerate  matrix.  Hence,  all  the  principal 
submatrices  of  M  are  also  nondegenerate.  Also,  every  subcomplementary  set 
of  column  vectors  is  linearly  independent  and  every  complementary  cone  has  a 
nonempty  interior. 

7.5  Let  {A  ^ . A.n-1^  any  su^comPlementary  set  of  column 

vectors.  We  will  now  show  that  the  hyperplane  generated  by  this  sub¬ 
complementary  set  of  column  vectors  strictly  separates  the  points 
representing  the  left  out  complementary  pair  of  column  vectors  I  and 
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Suppose  not.  Then  the  interiors  of  the  complementary  cones 
Pos{A^,  A>n_^,I.n>  and  Pos{A  A^n_^,-M  }  have  a  nonempty 

intersection. 


/ 


Then  we  pick  a  pair  of  points  q  ,  q  ,  both  nondegenerate  with  respect  to  M  , 

A 

sufficiently  close  together,  such  that  q  is  in  the  intersection  of 

Pos{A  . . A  ., I  }  and  Post  A  , ,  ...,  A  .  ,-M  }  and  q  lies 

•x  •n—x  *11  -n-l  *n 

A 

outside  both  these  cones  and  is  strictly  separated  from  q  by  the  hyperplane 

through  Pos{A  . . A  ,  }  .  We  can  choose 

• i  • n-i 

.  n-1 

q  -  I  A. A  .  +  ol 
J-l  J 


where  A.,  ...»  A  a  are  all  >  0  ,  a  sufficiently  small  A. . A  ,a 

x  n-x  l  n 

a  u“  X 

are  such  that  q  is  nondegenerate  with  respect  to  M  and  £  A  A  does 

J-i  J 

not  lie  in  any  subspace  generated  by  (n  -  1)  or  less  column  vectors  of 


(1,-M)  excepting  Pos{A^,  ...,  A^n_^}  .  By  the  nondegeneracy  of  q  , 

A 

q  +  AT  is  also  nondegenerate  with  respect  to  M  for  all  but  a  finite 
p  n 

number  of  values  of  A^  .  So  we  could  pick  an  >  0  sufficiently  small 


such  rhac 
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q  E  l  XjA.j  +  6I-n 

J-l  J 

„  nondegenerate  with  respect  to  M  for  all  6  satisfying  '  *  0  , 

<  B  <  a  .  So  if  we  had  chosen  our  original  a  so  small  that 

o  <  a  <  a  ,  and  q  is  in  the  interior  of  both  ?os{ktV  ••  A.n-l,I-n 
and  PostA^,  •••»  A.n-1,-M.n^  *  then 

-  n-1 

“  ‘  )h  V-J  ‘  0l-n 

ts  outside  both  these  complementary  cones  and  ",  does  not  lie  in  any 
complementary  cone  in  which  0  does  not  lie.  Thus,  the  number  of  solutions 
t„  (1)  when  ,  .  ;  U  strictly  less  (by  two,  than  the  n-ber  when  q  •  1  . 

leading  to  a  contradiction. 

By  e  similar  argument,  we  verify  that  the  hyperplane  through  any 
subcomplementary  set  of  colman  vectors  strictly  separates  the  points 
representing  the  left  out  complementary  pair  of  column  vectors. 
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7.6  We  will  now  show  that  the  principal  subproblem  of  (1)  in  the 

variables  (w. ,  w  ,  ;  z, ,  ...»  z  ,)  satisfies  a  similar  seperation 

property.  The  column  vector  corresponding  to  w^  in  this  subproblem  is 

the  j-th  column  vector  of  the  unit  matrix  of  order  (n-1)  ,  which  we  denote 

by  J  ,  and  the  column  vector  corresponding  to  z^  in  this  subproblem  is 

T 

-(m  ,m»  ,  . ...  m  .  .)  which  we  denote  by  -m  .  We  note  that  the 

lj  n-l,J  ’j 

column  vectors  in  the  subproblem  are  obtained  by  deleting  the  last  component 
from  the  column  vectors  in  the  original  problem. 

Let  {a  ...»  a  ^_^,a  ...»  a  ^  be  any  subcomplementary  set 

of  column  vectors  in  the  subproblem.  We  want  to  show  that  the  hyperplane 
in  Rn  *  through  these  column  vectors  strictly  separates  J  ^  and  -m^  • 
Let  A  ^  be  the  column  vector  corresponding  to  a  f,r«l,  •••»  i-l» 
i+1,  ....  n-1  ,  in  the  original  problem.  Then  (A  ...,  A  ^  ^,A 
....  is  a  subcomplementary  set  in  the  original  problem.  By 

7.5  the  hyperplane  in  Rn  through  these  column  vectors  strictly  separates 
and  -M^  .  Suppose  this  hyperplane  is 


DX  -  0 


where  D  -  (d^ . d^)  and  X  t  R  .  Then 


DA  .  -  0 

•  1 


DA.i-l  "  0 

DI  .  >0  and  D(-M  .)  <  0 

•  i  •  i 

BA.i+i  ■  0 


DA  .  -  0 
•  n-1 

DI  -  0  . 
•n 
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Now  DI  =  d  *=  0  .  Let  d  =  (d,  ,  ...»  d  , )  .  Because  d  =  0  ,  the 
•n  n  1  n-1  n 

above  equations  imply  that 


Ja.i-i  •  0 

dJ)  .  >0  and  d(-m  . )  <  0 

•  i 

da-i+i  '  0 


da  ,  *  0  . 
•n-1 


Let  x  ■  (x, . x  . )  e  Rn  *  .  Thus  dx  *  0  is  the  hyperplane  in  Rn  ^ 

1  n-l 

through  the  subcomplementary  set  {a^,  •••*  a.i_i»a.i+i»  •••»  a.n-l^ 
of  the  subproblem  and  it  strictly  separates  «J  and  -m  . 

Hence  the  subproblem  also  satisfies  a  similar  separation  property. 

By  a  similar  argument  we  can  verify  that  every  principal  subproblem  of  (1) 
of  order  (n-1)  satisfies  the  separation  property. 


7.7  Induction  Hypothesis: 

For  any  complementarity  problem  of  order  r  <  n-1  ,  with  column  vectors 
(I,-N)  ,  if  N  is  nondegenerate  and  if  the  hyperplane  through  every  sub¬ 
complementary  set  of  column  vectors  strictly  separates  the  points  representing 
the  left  out  complementary  pair  of  column  vectors  in  the  problem,  then  N 
is  a  P-matrix. 


7.8  The  induction  hypothesis  is  easily  verified  for  the  case  r  =  1  . 
By  nondegeneracy  N  ■  (ra^)  t  0  •  Since  r  ■  1  ,  the  sub complementary  set 
is  the  null  set  and  hence  the  hyperplane  through  the  subcomplementary  set 
is  the  singleton  consisting  of  the  origin  itself.  Since  this  separates 


[ 


the  points  on  R1  representing  1  and  -m^  *  we  should  have  -ra^  <  0  •  . 
So  N  =  (m.^)  >  0  and  hence  is  a  P-matrix  in  this  case. 

7.9  Hence  by  the  induction  hypothesis  and  7.6  every  principal  submatrix 
of  M  or  order  n-1  is  a  P-matrix.  Hence  all  principal  subdeterminants  of 
M  of  order  <  n-1  are  strictly  positive. 

Since  M  is  nondegenerate  by  7.4,  determinant  of  M  /  0  .  So  M  ^ 
exists.  So  the  constraints  (1)  can  be  written  as 


-1 

z  -  M  w  *  Q 

T 

z>0,w>0  ,  z  w  *  0 


where 


Q  -  -M_1q  . 


If  (1)  has  a  constant  number  of  solutions  for  every  q  e  Rn  ,  q  +  0  , 
then  (19)  has  a  constant  number  of  solutions  for  each  Q  e  Rn  ,  Q  /  0  . 

Hence  by  the  arguments  used  previously  all  principal  subdeterminants 
of  M  ^  or  order  (n-1)  or  less  are  strictly  positive.  Let  a  be  the 
value  of  the  principal  subdeterminant  of  M  *  obtained  by  striking  off 
the  first  row  and  column  from  M  ^  .  Then 


determinant  of  M 


But  a  >  0  ,  >  0  .  So  by  (20)  the  determinant  of  M  is  also  strictly 

positive.  So  all  principal  subdeterminants  of  M  are  strictly  positive. 
Hence  M  is  a  P-matrix. 

So  the  induction  hypothesis  7.7  holds  when  r  -  r  also.  It  has  been 
verified  for  r  -  1  in  7.8.  So  by  induction  it  holds  for  all  n  . 

Hence  by  7.5,  Theorem  7.2  is  true  for  all  n  . 


42 


8.  THE  ODD  NUMBER  THEOREM  FOR  NONNEGATIVE  Q-MATRICES; 

8.1  Here  we  show  that  if  M  is  a  nonnegative  Q-matrix,  then  the 
number  of  complementary  feasible  solutions  is  an  odd  number  whenever  q  is 
nondegenerate  with  respect  to  M  .  This  result  may  not  hold  if  M  is  not 
nonnegative . 

8. 2  Theorem: 

If  M  >  0  and  is  a  Q-matrix,  then  the  number  of  complementary  feasible 
solutions  Is  an  odd  number  for  any  q  nondegenerate  with  respect  to  M  . 

Ptoof : 

Proof  is  by  induction  on  n  . 

8.3  If  n  *  1  ,  then  M  *  (m^)  and  M  is  a  Q-matrix  if  and  only  if 
>  0  by  Theorem  5.2.  Here  q  ■  (q^  and  for  each  q  e  R1  there  is 

exactly  one  complementary  feasible  solution.  Hence  Theorem  8.2  is  true 
when  n  -  1  . 

8.4  Induction  Hypothesis: 

Suppose  Theorem  8.2  is  true  for  all  complementarity  problems  of  order 
(n-1)  or  less.  We  will  now  show  that  this  implies  that  Theorem  8.2  also 
holds  for  problems  of  order  n  . 

8.5  By  Corollary  5.4  all  principal  submatrices  of  M  are  also  Q-matrices. 

Consider  the  principal  subproblem  in  (w0,  ...,  w  ;  z.,  .  z  )  with 

i  n  l  n 

the  right  hand  constants  ■  Q.  .  If  Q,  is  nondegenerate  in  the  subproblem, 
then  it  has  an  odd  number  of  complementary  feasible  solutions  when  Q.  *  Q.  • 
by  the  induction  hypothesis  8.4. 
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-  /  ql  \ 

Let  q  *  [ - 1  where  q.  >  0  ,  be  nondegenerate  with  respect  to  M 

\  a  / 

Since  q^  >  0  and  M  >  0  ,  the  variable  z^  must  be  equal  to  zero  in 

and  solution  to  (1)  when  q  =  q  .  Thus  if  (w;z)  is  a  solution  to  (1) 

when  q  *  q  ,  z.  <=  0  and  hence  (w0,  w  ;  z z  )  is  a  complementary 

1  L  n  L  n 

feasible  solution  to  the  subproblem  when  Q  =  Q  . 

Also  if  (w*,  . ..,  w*  ;  z*,  . ..,  z*)  is  any  complementary  feasible 
solution  to  the  subproblera  when  Q.  “  Q.  »  define 


w*  ** 


>  0 


and  then  (vj,w*.  w*  ;  0,z*,  ....  z*)  is  a  complementary  feasible 

solution  to  the  original  problem  when  q  «  q  . 

Thus  uvery  complementary  feasible  solution  of  the  original  problem 
leads  to  a  complementary  feasible  solution  to  the  subproblem  and  vice  versa. 
Hence  boll  jiems  must  have  the  same  number  of  complementary  feasible 
solutions.  Hence  when  q  ■  q  ,  (1)  has  an  odd  number  of  solutions. 

By  a  similar  argument  we  conclude  that  the  original  problem  has  an 
odd  number  of  complementary  feasible  solutions  whenever  q  is  nondegenerate 
with  respect  to  M  and  at  least  one  component  in  the  vector  q  is  positive. 

It  only  remains  to  be  shown  that  the  same  result  holds  even  when 
q  <  0  . 

8.6  We  will  now  show  that  on  every  unbounded  edge  of  K(q)  lying  in 
the  almost  complementary  set  C^(q)  ,  both  the  variables  w^  and  z^ 

tend  to  4-®  while  zot  ....  z  remain  finite. 

2  n 

From  2.5,  K(q)  is  the  set  of  all  (w;z)  satisfying 


v  -  M  z  +  q  ,  i  -  1,  ....  n 

z  >  0  ,  w  >  0 

m  w  m 


(21) 


Since  M  is  a  Q-matrix,  m^  >  0  for  all  i  =  1,  ...»  n  by  Theorem  5.2. 
Consider  any  unbounded  edge  of  K(q)  .  If  all  the  variables 

z,  ,  ....  z  remain  finite  on  this  edge,  then  by  (21)  all  the  variables 

1  n 

w,  ,  ....  w  also  remain  finite  and  hence  the  edge  cannot  be  an  unbounded 
I  n 

edge.  Hence,  on  every  unbounded  edge  in  K(q)  at  least  one  of  the 

variables  z, ,  ....  z  must  tend  to  +°°  .  If  z.  tends  to  +«>  on  this 

edge,  then  from  (21)  and  the  facts  that  M  >  0  ,  m^  >  0  and  Is 

finite  and  fixed,  w^  must  also  tend  tc  +°°  along  this  edge.  Hence  if 

any  unbounded  edge  of  K(q)  lies  in  the  almost  complementary  set  C^(q)  , 

then  the  variables  z_,  ...»  z  should  all  remain  bounded  on  that  edge. 

2.  n 

Hence  z^  must  tend  to  +°°  on  that  edge  and  consequently  w^  also  tends 
to  +°°  on  that  edge. 

Thus  on  every  unbounded  edge  in  C^(q)  ,  the  variable  -*•  +®  . 

8.7  Suppose  q  is  nondegenerate  wJ  b  respect  to  M  .  Then  there  exists 

an  >  0  such  that  for  all  a  >  aQ  ,  the  point  q  -  al  ^  is  nondegenerate 
with  respect  to  M  .  Hence  the  entire  half-line 

{q  :  q  -  q  -  ol^  ,  a  >  aQ}  (22) 

lies  in  the  interior  of  a  set  of  complementary  cones.  We  now  show  that  the 

number  of  complementary  cones  in  which  this  half-line  lies  is  precisely  the 

number  of  unbounded  edges  in  C^(q)  .  Let 

F  ■  {(w;z)  :  (w;z)  -  (w*  +  ©w^jz1  +  0z^)  ,  0  >  0} 

be  an  unbounded  edge  in  C^(q)  .  Then 

(vj  +  0w*)(zj  +  0z*)  -  0  for  all  i  +  1 

2 

and  w^  >  0  by  8.6.  Hence 


for  all  0  >  0 
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/12  12  12  i  2  \ 

(v;z)  -  ^0,w2+0w2,  ....  wn+6wn  ;  z^+Bz^  ....  zn+ezn/ 

is  a  complementary  feasible  solution  for 

q  *  q  -  +  0w^jl  ±  for  6  >  0 

2 

and  since  w^  >  0  ,  as  0  varies  from  0  to  ® 

jq  :  q  »  q  -  (vj  +  Ow^I^  .  6  >  °} 

is  eventually  the  same  half-line  as  in  (22). 

Also  for  any  a  >  c»o  ,  q  -  al  ^  cannot  lie  in  any  complementary  cone 
which  has  Pos{I  as  a  generator.  For,  if  it  does,  there  exists  a 
subcomplementary  set  of  columns  {B  2 . B^}  suc^  that 

-  n 

q  -  al  ,  -  A  I  +  l  X  B 

i-2 

for  some  X, ,  ...,  X  >  0  .  Then  q  -  (a  +  X.)I  .  lies  in  the  subspace 
in*  x 

through  the  subcomplementary  set  {B>2,  •••»  B>n)  contradicting  the 
assumption  that  q  -  al^  is  nondegenerate  with  respect  to  M  for  all 
a  >  a 

o 

Hence  if  the  half-line  in  (22)  lies  in  some  complementary  cone,  say, 

Pos{A'ltA'2»  ...»  A>n)  then  A#1 . A>n  must  be  linearly  independent 

and  A  ,  »  -M  ,  .  Then  we  can  express  this  half-line  as 
•1*1 

f  !  ’  ■  Jl6!4'1  +  (a  "  ^  ’  0  ”  “°} 

1  2 

for  some  B  ,  B  >  0  .  Thus 


“*“*  *  “"WWUMWSVi' 
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q  =  al  +  l  6jA  +  (a  -  a  )  [  B*A  (23) 

i-1  i»l 

for  any  a  >  . 

r  r  12 

Suppose  (w  ;z  )  is  obtained  by  setting  «=  1  ,  and  the 

variable  associated  with  the  column  vector  A  ^  equal  to  B^  ,  i  *  1 . n 

and  all  the  other  variables  in  (w;z)  equal  to  zero,  for  r  ■=  1,2  .  Then 

(22)  implies  that 

F  *  {(w;z)  :  (w;z)  =  (w^  +  0w^  ;  z^  +  0z^)  ,  0  >  0} 
is  an  unbounded  edge  in  C^(q)  . 

Thus  every  unbounded  edge  in  C^(q)  gives  raise  to  a  complementary 
cone  in  the  interior  of  which  the  half-line  in  (22)  lies  and  vice  versa. 

Hence  the  number  of  unbounded  edges  in  C^(q)  is  equal  to  the  number  of 
complementary  feasible  solutions  for  q  -  al  ^  where  a  is  a  sufficiently 
large  number. 

*  *»  -  ip 

8.8  Thus  for  any  q^  such  that  q  -  (q^^,  ....  q^)  is  nondegenerate 

a 

with  respect  to  M  ,  the  number  of  unbounded  edges  in  C^(q)  is  a  constant. 
This  number  is  equal  to  the  number  of  complementary  cones  in  which  the  half¬ 
line  (22)  eventually  lies  as  aQ  is  made  large. 

8.9  By  the  nondegeneracy  of  q^  we  know  that  there  exists  a  Bq  such 

-  x 

that  for  all  8  >  Bq  ,  (B,q2»  •••»  9n)  is  nondegenerate  with  respect  to  M  . 

'  -  t 

Hence  we  can  always  pick  a  q*  >  0  such  that  q*  *  (q*,q2>  •••»  Tn)  is 
nondegenerate  with  respect  to  M  .  Since  q*  >  0  ,  the  number  of  complementary 
feasible  solutions  when  q  -  q*  is  an  odd  number.  Therefore  by  6.3(iii)  the 
number  of  unbounded  edges  in  C^(q*)  is  an  odd  number  and  hence  by  8.8  the 
number  of  unbounded  edges  in  C^(q)  is  an  odd  number.  By  6.3(iii),  the 
number  of  complementary  feasible  solutions  when  q  -  q  is  therefore  an  odd 


number. 
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Hence  under  the  induction  hypothesis,  Theorem  8.2  holds  for  the  original 
problem  of  order  n  .  By  8.3  and  by  induction,  Theorem  8.2  is  true  for  all  n  . 

8.10  Corollary: 

If  M  is  a  Q-matrix  and  if  there  exists  a  complementary  set  of  column 

vectors  {A  ^ . A^}  is  linearly  independent,  such  that  each  of 

the  remaining  vectors  B  ,  ...,  B  among  the  column  vectors  of  (I,-M) 
satisfies 


Bj  c  Pos{-A#^,  ...»  -A>n)  for  all  J  ■  1,  ...,  n 

then  the  number  of  complementary  feasible  solutions  is  an  odd  number  for  all 
q  which  are  nondegenerate  with  respect  to  M  . 

Proof : 


Transform  the  column  vectors  A  ,,  ...,  A  into  the  column  vectors  of 

•  1*  ‘n 

the  unit  matrix  by  making  the  necessary  principal  pivots.  Then  Corollary 
8.10  follows  from  Theorem  8.2  and  2.21. 

8.11  In  the  special  case  when  n  ■  2,  the  restriction  that  M  >  0  can 
be  removed  from  the  hypothesis  of  Theorem  8.2.  This  is  discussed  below. 


8.12  Theorem : 

If  n  *  2  and  M  is  a  Q-matrix  then  the  number  of  complementary 
feasible  solutions  is  an  odd  number  whenever  q  is  nondegenerate  with  respect 
to  M  . 


Proof : 


8.13  Case  1:  If  Posf-M^.-M^}  is  a  subset  of  the  nonpositive  orthant 
2 

of  R  ,  Theorem  8.12  follows  from  Theorem  8.2. 

8.14  Case  2:  If  Pos{-I  .  .-Io)  C  Pos{-M  ,  ,-M  „}  the  hypothesis  that 

— — —  .X  *■ 

M  is  a  Q-matrix  implies  that  -M  ^  and  -M  2  are  contained  one  each  in 
Pos{I#1,-I<2}  and  Pos{-It^,I  2>  respectively. 


1 


h5 Scoria  ft 


We  verify  that  in  this  case  the  number  of  complementary  feasible  solutions 
is  1  or  3  for  every  q  nondegenerate  with  respect  to  M  . 


8.15  Case  3;  Since  M  is  a  Q-matrix,  the  only  other  possibility  is 

that  exactly  one  of  -M  or  -M  is  contained  in  the  interior  of 

•  1  •  2 

Pos{-I  ,-I  }  .  Suppose  it  is  -M  .  Then  the  hypothesis  that  M  is  a 

•  1*2  *1 

Q-matrix  implies  that  either  -M^  e  Pos^1  .i»M.i^  or  _M-2  e  Pos^ .1,_1 .2^  ’ 


In  either  case  we  verify  that  the  number  of  complementary  feasible  solutions 
is  either  1  or  3  for  all  q  nondegenerate  with  respect  to  M  . 


KMUrVMBIM 


hwawswfiti  ik«* 
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8.16  Corollary: 

If  n  =  2  ,  there  exists  a  q  nondegenerate  with  respect  to  M  ,  for 
which  the  number  of  complementary  feasible  solutions  is  at  most  one. 

8.17  Note:  When  n  >  3  ,  Theorem  8.2  is  not  necessarily  true  if 
M  |  0  ,  and  Corollary  8.16  may  not  be  true. 

As  an  example  consider 


It  can  be  shown  that  this  is  a  Q-matrix  by  verifying  that  the  union  of  all 

3 

the  8  complementary  cones  is  R  .  Also  M  is  a  nondegenerate  matrix.  We 

T 

verify  that  q  -  (1,1,1)  is  nondegenerate  with  respect  to  M  .  When  q  ■  q 
there  are  four  distinct  complementary  feasible  solutions,  because  q  lies 
in  each  of  the  complementary  cones  Pos{I  *  P°s(-M  ^,1  * 

Post  I  .  ,-M  _,I  «}  and  Post  I  .  ,1  _,-M  .}  and  in  none  of  the  others. 

By  Theorem  6.2,  the  number  of  complementary  feasible  solutions  is  an 
even  number  for  all  q  nondegenerate  with  respect  to  M  ,  and  since  M  is 
a  Q-matrix,  this  number  must  be  >  2  . 

This  shows  that  the  converse  of  Corollary  6.9  is  not  necessarily  true 
unless  M  >  0  . 

8.18  Note:  When  n  >  3  ,  the  number  of  complementary  feasible  solutions 
can  be  2  for  all  q  c  Rn  .  The  example  in  8.17  shows  this.  Thus  when 
n  >  3  the  complementary  cover  can  span  the  whole  space  more  than  twice 


around. 
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